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We present an example of a gauge-invariant Lagrangian that contains four derivatives and de- 
scribes one massive, non- ghostlike, degree of freedom. 
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Motivation. Gauge invariance is often assumed to 
imply the existence of massless degrees of freedom, 
although there are famous counter-examples, ranging 
from 1, 2] to 3, 4]. Moreover, it is often assumed that 
theories with higher derivatives contain ghosts or other 
pathologies. Recently, [f| has observed that the model, 
proposed in Q and therein called "new massive grav- 
ity," does provide a counter-example to both these com- 
monly accepted no-goes. Alas, the model |6| is defined 
only in three-dimensional space-times. 

In the present paper we will discuss a model that vi- 
olates both aforementioned assumptions in any number 
of dimensions. As we will see, the model shares many of 
the properties with the new massive gravity of [6] , but 
without the limitation of three space-time dimensions 
(though we will mostly focus, for the sake of argument, 
on four dimensions). 

The key ingredient of this model is, in d space-time 
dimensions, a id — l)-form field potential Ap 1 p 2 ...p d _ 1 . 
As for electromagnetism, which is described by a 1-form 
potential, we require gauge invariance under the trans- 



formation A 



Pip.2---p.d-l 



A 



MlM2---Md-l 



- i s an arbitrary (d — 2)-form. Thus, the 



is gauge 



where B t 

field strength F^i^-tid-i = 5 r„*4 
invariant. 

Specializing to four space-time dimensions, it is usu- 
ally assumed, by analogy with electromagnetism, that 
the Lagrangian for a 4-form field strength should be 



1 



-pvpX 



T, 



(1) 



Had this Lagrangian - or, more generally, any polyno- 
mial in Tp Vp x - been chosen, the field Ap Vp would not 
have any local dynamics. We observe that the tenso- 
rial structure of the field strength imposes Tp Up x(x a ) — 
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yf— 1)£pvp\ $(x a ), where <&(x a ) is a scalar field and 
£pvp\ the Levi-Civita symbol with £0123 = 1- The equa- 
tions of motion derived from ([T]), V^T^upX = 0, imply 
<9 M $ = 0, i.e., $ is constant. 

We will show that adding higher derivatives term, 
V Q J-" Ml/p A, in the Lagrangian allows the (dual) scalar 
field $ to propagate. Moreover, for an appropriate 
choice of parameters, the propagating degree of freedom 
has a well-behaved kinetic term and a positive-definite 
potential. 

The Lagrangian. The most general Lagrangian, 
quadratic in T and with at most four derivatives 
and "dimension six" operators, is (we will assume the 
"mostly plus" convention for the signature of the met- 
ric) 



48 ^ pA+ 48M 2 



Cl 



V a ^ pX VaFp»pX + 



C-2 



12 M 2 



VaT^^V^T^px + terms in RTT 



(2) 



where RTT indicates any term built by contracting two 
field strengths with a Riemann tensor or its contrac- 
tions. We will neglect these terms to ease the discus- 
sion, but their role will be discussed at the end of this 
note. In ([2]), the coefficients cq, ci, C2 are dimensionless 
numbers, and M is a mass scale. A term proportional 
to V a J r/3l/pA V pT avp x can also be added to (HJ, however 
it is identical to the term proportional to C2 after an 
integration by parts and so it will be neglected. 

It is convenient to rescale Ap Vp by a factor of M; from 
now on, we will work with the field Ap Vp = Ap Vp /M 
(and F^px = Tp Vp x/M). The equations of motion 
derived by varying the Lagrangian ([2]) with respect to 
A vp x can be rearranged to take the form 

V p [co M 2 F^x - (ci + c 2 ) \7 a \7 a Fp Up x] = , (3) 

where we have used the fact that Fp Vp x is totally 
antisymmetric in its four indices. Setting Fp Up x = 
\/—g tpupx 4>, the equation of motion becomes 



(ci + c 2 ) V Q V Q 4> - c M 2 4> = g , 



(4) 



2 



where q is an integration constant. Therefore, the La- 
grangian ([2]) describes one dynamical scalar degree of 
freedom of mass 



Co 



ci + c 2 



M 



(•5) 



and vacuum expectation value proportional to the inte- 
gration constant q. In other words, our system is equiv- 
alent to a landscape of massive scalars characterized by 
different values of the integration constant q. 

The stress-energy tensor obtained by varying the La- 
grangian ([2]) with respect to reads, in terms of the 
scalar <b. 



T^v = (ci + c 2 ) 
or, by using (j4]), as 



^ -(ci +c 2 ) V^V^ 



-Paw 



(0) 



(7) 
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The requirement that the stress-energy tensor is pos- 
itive definite imposes the conditions c\ + c 2 > (equiv- 
alent to the no-ghost condition for the scalar kinetic 
term) and c > (equivalent for the mass of 4> to be 
non-tachyonic) . 

We note that the requirement of non-tachyonic mass 
term corresponds to the "wrong" sign for the F 2 term. 
This is precisely the same phenomenon observed in the 
new massive gravity model of @, where the role of F 2 
is played by the Ricci scalar. Let us also note, however, 
that the addition in the Lagrangian of higher powers of 
the form field, e.g., F a ^ s F p » pX F aPpX F p ^ 5 , can stabi- 
lize the vacuum of the theory at non- vanishing values 
of the form field even if cq < 0. 

The properties of our system become more apparent 
in the first order formalism, where the potential A pvp 
and the field strength F pup \ are treated as independent 
variables, and the condition F^ vp \ = d[ p A up x\ is im- 
posed by introducing a Lagrange multiplier q(x a ). The 
Lagrangian then reads 



co M 2 



F^ pX F pupX + -L\7 a F^ px V a F^ pX - 
4o 



(8) 



Since F^ vpx is now a fundamental field, we can re- 
place it by \/—ge pV p\ 4>{x a ) directly in the Lagrangian: 



ci + c 2 



co 
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AivpX a A 



(9) 



The equation for A p,vp , d\q = 0, simply implies the 
Lagrange multiplier q to be a constant. 

Despite this system being very similar to the one de- 
scribed in 0, IH, it is not the same: while in 0, Q the 
different vacua have all the same cosmological constant, 
in our case they do not, since the vacuum energy is pro- 
portional to q 2 . Further differences emerge when we 
couple our form field to matter. 

Coupled to matter. Given the "wrong" sign of the 
term F 2 , an interesting question to answer is whether 
this model is stable once coupled to matter, p-form 
fields couple naturally to (p— l)-branes. In four dimen- 
sions, a 2-brane sources a 3-form field 



d 3 C 



: dz p dz v dzP 

: > pvp d^WW e 



abc 



(10) 



where e is the charge and a the tension of the mem- 
brane; we denote by d the brane world-volume. Note 
that e is proportional by a factor M to the original 
brane charge for the field A pvp . This coupling could 
destabilize the vacuum through brane nucleations, de- 
scribed by the following equations of motion (neglecting 
the role of gravity) 



v v 8 (VM - —(h^ w P x dZp dZp dzx r abc 



(11) 



U<j>- fi 2 (j)- 



Q 



c\ + c 2 



0, 



(12) 
(13) 



where we have already substituted the 4-form curvature 
of A pvp with F^px — e^ p \ <j) an( l is the effective 
mass found in eq. ([5]). 

The solution for the equation of motion for the 3- 
brane (|11[) is that of a hyperboloidal bubble of con- 
stant "radius" r$. From (jT2"j) . we notice that the La- 
grange multiplier q is constant away from the bubble 
and changes by e, the charge of the 2-brane, when cross- 
ing the bubble, i.e., |Ag| = e. 

Since we are interested in instanton solutions for the 
nuclcation of bubbles with different value of q, we will 
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consider the Euclidean action Se of ^ with coupling 
(fTOf on the solution of the equations of motion (fTTj) . 
CU) and (USD 



5 E 



d 4 x 



+ 



(14) 



where is the canonically normalized scalar field <j> = 
c0 with c = \/ci + C2, and where we have taken into 
account the appropriate boundary terms in the La- 
grangian. The function q(r) is given by q(r) = q < 0(ro — 
r) + q > 6{r — ro), where the step function 8(r) is equal 
to 1 when its argument is positive and equal to zero 
when it is not, and |g> — q< \ = e. By using the spheri- 
cal symmetry of the problem, we can simplify (|14[) even 
further 



dr r A 



(d r & 2 + /*V + 2 



g(0 | 
c 



2 n 2 a r% 
(15) 



and the equation of motion for <f> becomes 



„</>' 2 7 i( r ) n 

o n <p — = . 

r c 



(16) 



Regularity for both r — > + and r — > oo constrains the 
solution to be 



In the first limit, the difference in actions reduces to 



A<S E - 2 tt 2 a r% 



7r 2 ?o <7>ee 



2^ 2 



(21) 



where we have introduced the symbol e = ±1, q> — g< = 
£ e, to discriminate between a brane and an anti-brane 
nucleation. This expression for the change in action can 
be easily found by direct calculation, observing that, 
since the bubble is much smaller than the scalar Comp- 
ton wavelength, (f> stays essentially constant at the value 
it assumes outside the bubble, i.e., <j> = —q > /(fj 2 c). 

It is clear that the action (f2"Tj) can always be mini- 
mized, with an appropriate choice of the sign e, at 
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(22) 
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In the opposite regime /iro ^ 1, the action reads 



2 a- r£ 



tt 2 g> £ e — e z 
4/i 2 c 2 



2/r 3 52 r °' 



(24) 

where the second term corresponds to the change in 
bulk energy induced by the bubble nucleation, whereas 
the third term represents a correction to the effective 
membrane tension induced by the gradient in the field 



(r) 




hXt±r} 
fi r 



for r < ro 
for r > tq 



(17) 



where I n (/ir) and K„(/ir) are Bessel functions. The 
constants A and B are found by requiring continuity 
of the function and its derivative across the bubble at 
r = r : 



A = 
B = 



1< 



— rlK 2 {fir ) 



g< ~ g> 2, I \ 



(18) 
(19) 



The difference between the action for the solution of 
constant q > and the action for q(r) is 



2 TT O fn 



7T 2 r% ql 



1% 



4/x 2 

K 2r r\ (q< - q>Y 



h(»r )K 2 ((ir ). (20) 



This expression is not very illuminating; it is conve- 
nient to consider the limits /iro "C 1 and /iro ^ 1. 



The bubble radius is found by extremizing the action 
with respect to r . Inspection of eq. (I2U1) , and of its ex- 
plicit limits (12H and (|24| , shows that it is always possi- 
ble to find a ro > that minimizes the effective action. 
This implies that, unlike the usual Brown- Teitelboim 
model 0, [l(| , this system does not have any absolutely 
stable vacuum. The reason for such an instability lies in 
the "wrong" sign for the form kinetic term, F 2 . Never- 
theless, it is easy to find parameters for which the decay 
rate, proportional to e~ SE , is extremely small and the 
system can be considered stable for all practical pur- 
poses. 

In the decoupling limit for the field <fi, M — > oo with 
q/M = q and e/M = e fixed, the action describing 
membrane nucleation reads 



4 c 



^ (2 q> £ e - e 2 



(25) 



which gives the Brown- Teitelboim result [tj H3] when 
c = -1. 

By solving numerically the equation d ro ASE = 0, 
we find that the probability to nucleate a bubble of 
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q < = q > — e is always smaller than the one correspond- 
ing to the Brown- Teitelboim limit, M — ¥ oo. This can 
be understood by considering that, when the field cf) is 
dynamical and massive, a part of the energy gained by 
nucleating a bubble goes into exciting </>-modes and the 
smaller the mass, the more likely these modes are to be 
excited. 

Equivalently, it can be noted that <f> generates an ef- 
fective attractive force between branes. For two parallel 
branes with charges e and e' placed at a distance d, the 
field <p mediates a force per unit surface 



the mass term for the vector field realized via the 
Schwinger mechanism effectively generates a kinetic 
term for the Lagrange multiplier Q, making it dynami- 
cal. The limit e — > decouples Q from the other fields; 
the condition d^Q = is recovered and Q ceases to 
propagate. 

Before canonically normalizing the fields, we notice 
that, for ci + C2 > 0, the sign of the kinetic term for 
both $ and Q is that of a healthy scalar field, that is 
not a ghost. 

Let us substitute (f> = V c i + c 2 $/M and q = \/nQ/e 



F 

s 



2 /i c 



— exp [-fi d] 



(26) 



that is, two branes of same charge repel each other, 
two of opposite charge attract each other. To defy this 
force the branes will have to be nucleated at radius tq 
greater than when such a force is absent, i.e., /i — > oo, 
the Brown- Teitelboim case. A larger bubble will be less 
likely to nucleate than a smaller one. 

Schwinger model with a dynamical two-form. The 
two dimensional version of our scenario is especially in- 
teresting. In fact, in two dimensions, it is known, Q, 
that once a massless fermion - the analog of a light 
brane in higher dimensions - of charge e is coupled to 
a vector field, a gauge-invariant mass term is generated 



^mass Afi 

Z 7T 



/!// 



-Q2~ 



(27) 



This mechanism works also when higher derivatives 
operators are considered - as it can be easily checked 
- therefore we simply add this gauge-invariant mass to 
the two dimensional version of the Lagrangian @. The 
longitudinal mode of A^ drops out, since the entire La- 
grangian is gauge invariant and it depends only on the 



transverse mode Aj, 



i- 4 -r^v 



AM' 



■(V Q J^) 2 



— (v a j- QAl ) 2 - 



2M 



„-A T u A T » 



— e M "C T 
2 fc K-'nv 



2V£), (28) 



where we have added a Lagrange multiplier Q to treat 
the vector field and the curvature independently; Tp, v is 
proportional to the two dimensional Levi-Civita symbol 
{T^v = e^v^) and A^ can be integrated out in favor of 
Q, Aj_ — —(ir/e 2 ) e iia d a Q. These substitutions lead to 
the following effective Lagrangian in terms of the scalar 
fields $ and Q, 



C 



off 



ci + c 2 
2M 2 



(d a <5>f 



2 



7T 

2~? 



(d a Q) 2 ; 
(29) 



e fi - 



+ 



(30) 



where fj, is the effective mass defined in (JSJ) - The model 
does therefore contain two propagating degrees of free- 
dom of masses: 



2 



1± 



4e 2 

7TC /i 2 



(31) 



Because the coefficient cq and the effective mass n 2 are 
positive definite, the square root is always greater than 
1, hence one scalar mode is always tachyonic. This 
tachyonic instability is the perturbative counterpart of 
the non-perturbative instability for brane nucleation de- 
scribed above. 

In the limit /i 2 3> e 2 /co, the two masses are approxi- 
mately fj, 2 and — e 2 /(7rco) and the tachyonic instability 
has a characteristic time much longer than timescale 
jjT . In the opposite regime [j? -C e 2 /c , the mass 

~), signal- 



eigenvalues have equal magnitude ±2 e /i/-^/ 
ing a fast instability. 

Non-minimal coupling to gravity. Let us then dis- 
cuss the possibility that the Lagrangian © contains 
couplings of T^px to the Riemann tensor. The formal- 
ism, that led to the equivalent expression ©, shows 
that such extra terms would sum up to a coupling pro- 
portional to <f> 2 R, with R, the Ricci scalar. Therefore, 
the effective Lagrangian for the canonically normalized 
field d> will read 



Mp 



(32) 

where the value of the parameter £ depends on the 
coefficients of the RJ- T terms in the Lagrangian ([2]). 
In ([32"T) we have also included the Einstein-Hilbert term 
for gravity, while we have already set q to be a constant, 
hence the term proportional to q d^A vp \ is neglected. 
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Direct couplings of J-^ up \ to the Riemann tensor, 
therefore, lead to a non-minimally coupled scalar field. 
The behavior of 4> is more transparent if we work in the 
Einstein frame. In order to do so, we redefine the met- 
ric = Q 2 g^, n 2 = 1 - £ ^ 2 /M 
Lagrangian is 



p. The transformed 



M 2 



R — 



n 2 + 6e 4> 2 /m 2 



v a 4>v a 4>~u{4>) 



2 " 2 ft 4 

(33) 

where U(<f>) = (fi 2 (j) 2 /2 + q (j)/d) /fl 4 . The formulae 
simplify for £ = 1/6, in which case the canonically nor- 
malized scalar reads 4> c = \/6 Mp arctanh f^/v6Mp), 
and has potential 



\/6Mj 



2 — / — / / 2 2 — / 

-) =cosh (f> c smh(/) c (3 /i Mpsinh^> c - 



- y/QMp cosh< 
c 



(34) 



The expression Q34[) shows that, for non- minimal cou- 
pling to gravity, the dependence of the system on the 
integration constant q becomes nontrivial. In this case, 
for instance, when q is large enough the potential be- 
comes unbounded from below. 

Discussion and conclusions. We have shown that 
the Lagrangian ([2]) is, despite its higher derivatives, a 
well-behaved Lagrangian: neither ghosts nor tachyons 
are propagating degrees of freedom. 

When we couple the form field to matter, instabilities 
do emerge, but, since they are non-perturbative, their 
characteristic time can be made safely long. Given the 
similarities between our system and the new massive 
gravity of [6j , it would be interesting to inquire whether 
analogous instabilities appear in that model when grav- 
ity is coupled to matter. As we noted above, the absence 



of an absolutely stable minimum lies in the "wrong" 
sign of the would-be kinetic term for the cZ-form, i.e., 
F 2 ; similarly in the model of @, the would-be kinetic 
term for the gravitational degrees of freedom, i.e., R, is 
chosen to have the "wrong" sign in order for the mass of 
the propagating degree of freedom to be non-tachyonic, 
hence we can speculate that a non-perturbative insta- 
bility similar to the one we found for the present model 
could appear when the new massive gravity is coupled 
to matter. 



As in our case, [ll| has also recently obtained a 
healthy behavior for form fields despite a higher deriva- 
tive Lagrangian. Our mechanism is however different, 
as one can see by observing that 



11| considers form 



fields whose rank is different from the one we are con- 
sidering. 

Let us conclude with a speculation. The goal of this 
paper is to show an example of a higher derivative gauge 
theory that does not have any unhealthy behavior and 
that propagates massive degrees of freedom. Our goal 
is not to embed this theory in a UV-complete model; 
it would be interesting, however, to study whether the 
Lagrangian ^ can emerge as a limit of some known the- 
ories, e.g., as an effective low-energy description arising 
from strong dynamics at scale M. Within this frame- 
work, the model we presented might be a generalization 
of the models describing a condensation of topological 
defects, as studied in \12\ . 
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